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We show how the dynamical modulation of the qubit-field coupling strength in a circuit quantum electrody¬ 
namics architecture mimics the motion of the qubit at relativistic speeds. This allows us to propose a realistic 
experiment to detect microwave photons coming from simulated acceleration radiation. Moreover, by com¬ 
bining this technique with the dynamical Casimir physics, we enhance the toolbox for studying relativistic 
phenomena in quantum field theory with superconducting circuits. 


Circuit quantum electrodynamics (cQED) m has swiftly 
become one of the leading quantum technologies for testing 
fundamentals of quantum mechanics, and for the implemen¬ 
tation of quantum information tasks. The possibility of fast 
tuning characteristic parameters of the interaction between 
superconducting qubits and the electromagnetic field allows 
for the experimental investigation of quantum optical systems 
in previously inaccessible regimes, thus enabling the explo¬ 
ration of a wide variety of new physical phenomena urn 
In particular, cQED is a natural arena to implement quantum 
field theory (QFT) concepts 0. Indeed, the physics associ¬ 
ated to the dynamical Casimir effect (DCE) —the generation 
of photons out of the quantum vacuum through the motion 
of boundary conditions—has been implemented in a super¬ 
conducting circuit platform. This experiment opened a new 
avenue of research on the analysis of the role of relativistic 
motion in quantum information setups 00. In this sense, 
quantum technologies are rapidly progressing from tabletop 
experiments to long-range networks ITOl and may enter, in the 
near future, into the space-based realm eh. In this scenario, 
relativistic effects are expected to become relevant 03 . 

The Unruh effect Q3], a paradigmatic quantum field the¬ 
ory prediction, consists in the excitation of ground state atoms 
due to their acceleration through vacuum, followed by spon¬ 
taneous emission. Such a phenomenon has never been ob¬ 
served, requiring unachievable accelerations in order to gen¬ 
erate detectable signals. This effect can be enhanced by many 
orders of magnitude when two-level quantum systems are ac¬ 
celerated through a single-mode high-Q cavity lfT4l . The rela¬ 
tion between this cavity-enhanced acceleration radiation and 
the original formulation of the Unruh effect in free space has 
produced a thought-provoking debate GSHHl. 

In this Letter, we show how the ultrafast variation of the 
coupling strength between a superconducting qubit and a 
transmission line resonator mimics the motion of the qubit at 
relativistic speeds. This analogy paves the way for novel ex¬ 
perimental implementations of QFT effects associated to rel¬ 
ativistic motion. Here, we propose the measurement of the 
radiation produced by an artificial atom performing simulated 
acceleration through a superconducting cavity. Our technique 


provides a novel tool for the quantum simulation of QFT in 
curved space times with current circuit QED technology, in 
addition to the observation of the DCE f4j and the possibility 
of generating effective spacetime metrics sin. 

We consider the system depicted in Fig. [T] consisting of 
a superconducting qubit interacting with the electromagnetic 
held confined within a transmission line resonator (TLR). 
This model can be described by the Hamiltonian 

H = tva)a + yu z + 'Hi(x q ), (1) 

where uj q is the qubit energy spacing, a z and a x are the usual 
Pauli operators acting on the qubit Hilbert space, and we con¬ 
sider h = 1. We assume that the system dynamics effectively 
involves a single resonator mode, described by annihilation 
and creation operators a and a', respectively, of frequency 
to = vk and wave vectors k = it/L. Here, L is the resonator 
length and v is the group velocity, which is established by 
the resonator electrical properties. The interaction Hamilto¬ 
nian between superconducting qubits and a transmission line 
resonator (TLR) is mediated either by the electrical poten¬ 
tial V, for capacitive coupling like in the case of transmon 
qubits, or by the current I for inductive coupling in the case 
of flux qubits. For the resonator stationary modes, both physi¬ 
cal quantities have a periodic space modulation. Without loss 



FIG. 1: (Color online) Sketch of a transmission line resonator cou¬ 
pled to a superconducting qubit, where the coupling strength depends 
on the position of the latter. The red dashed line follows the voltage 
profile of the second resonant mode of a A/2 resonator. 
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of generality, we can write the interaction Hamiltonian as 

Ui(x q ) = gcos(kXq)<j x (a) + a), (2) 

where g is the coupling strength and x q the qubit position fl9l . 
We point out that in circuit QED designs, it is possible to 
modulate g in time, by means of a controllable magnetic flux 
<t>x = / 0o flowing through the loop of a superconduct¬ 
ing quantum interference device (SQUID). We introduced the 
frustration parameter / and the magnetic flux quantum Dq. 
Due to the non-linear inductance of the Josephson junctions, 
such tunable coupling has usually a periodic dependence on /, 
g = go cos (/). We consider now a fixed qubit position, cho¬ 
sen in order to maximise the coupling strength. This is done 
by placing the qubit position in a node(antinode) of the cur¬ 
rent I for a charge(flux) qubit, corresponding to x q = 0 in 
Eq. 0 Accordingly, the interaction Hamiltonian has the fol¬ 
lowing dependence on /, 

'kLi(f) = 9 o cos (f)a x {a^ + a). (3) 

Choosing the profile for the magnetic flux 

f = k Xq, ( 4 ) 

the Hamiltonians in Eqs. 0 and 0 are equivalent. Therefore, 
the modulation of the effective coupling constant mimics the 
motion of the qubit x q (t ) inside the TLR. Later on, we will de¬ 
tail how this effective motion can produce qubit and field ex¬ 
citations in the relativistic regime. In the following, we show 
that fast-tuning of the coupling strength can be achieved in the 
framework of superconducting circuits. 

There are different ways of implementing a qubit-resonator 
system interacting through the tunable Hamiltonian of Eq. 0, 
using either phase |20]| . flux ETl l22ll or transmons |[23l l24l 
qubits. The configuration here considered, shown in Fig.0 al¬ 
lows independent tuning of the qubit-field interaction g and of 
the qubit energy spacing uj. Hence, it provides a versatile so¬ 
lution for simulating a broad range of parameter regimes with 
a single sample. As shown in Fig. 0a), the present qubit con¬ 
sists of two capacitively shunted SQUID loops, which form 
three superconducting islands. This design is well described 
with a model composed of two transmons, i.e. anharmonic os¬ 
cillators, capacitively coupled with the TLR and between each 
other. The qubit will be formed by the two lowest collective 
eigenstates of such system, here labeled with \Eq) and \Ei). 
As shown in Fig. 0b), when the two anharmonic oscillators 
are detuned, the qubit logical states are given by the vacuum 
|0 )l = I Eo) = |0) a |0) 6 and |1) L = \E X ) = |l) o |0) 6 , where, 
without loss of generality, we assume that u a < u>b- We de¬ 
note with |0)j and |l)j the ground and first excited eigenstate 
of the ?'-th transmon, and with to, the energy spacing between 
such states. On the other hand, when the two anharmonic os¬ 
cillators are nearly degenerate, the first collective excited level 
will be given by |1) L = |£i) = (|l) a |0) 6 - |0) a |l) 6 ) /V2- 
Such a state has no dipole coupling with the TLR, so in the 
latter configuration the cavity-qubit interaction is strongly in¬ 
hibited l23ll . 



FIG. 2: (Color online) (a) Sketch of a modified transmon qubit, see 
Ref. (231. Different colors identify independent superconducting is¬ 
lands capacitively coupled between each other, (b) Scheme of the 
energy levels. The qubit is composed of two coupled transmons, the 
logical states are given by the two lower collective eigenstates (blue 
and red). The third collective energy level (green) can be used for 
reading the qubit state, (c) Sketch of a superconducting chip con¬ 
taining a TLR interacting with a tunable coupling qubit. 


Controlling the external magnetic flux flowing through the 
two SQUID loops allows to independently tune the level spac¬ 
ing of the two transmons. This enables us to swap between the 
two collective energy configurations and, at the same time, to 
control the logical qubit transition frequency. Using this con¬ 
figuration, tuning of the cavity-qubit coupling strength from 
hundreds of KHz to tens of MHz has been experimentally 
proven j25l . We also point out that, in the nearly degenerate 
case (ui a ss uib), the second excited collective level is given by 
|2) = I E 2 ) = (|l)a|0)b + |0) a |l)&) /y/2. This state experi¬ 
ences a strong dipole coupling with the cavity mode, hence it 
can be used for fast-reading the qubit state. Such a measure¬ 
ment scheme consists in driving the transition \E{) —> \E 2 ) 
and then measuring the shift of the cavity resonant frequency 
due to the qubit state. Notice that the transition to the third 
level is safely detuned by at least 1 GHz from the primary 
qubit transition l25l . 

On-the-fly tuning of the coupling strength over nanosecond 
time scales has already been demonstrated in circuit QED ar¬ 
chitectures ll26ll . Fast modulation of the cavity-qubit coupling 
strength has also been proposed for similar designs |27l . In 
our scheme, the effective simulation of the qubit position is 
achieved through the modulation of the magnetic flux thread¬ 
ing SQUIDs, that is, similarly to the motion of the boundary 
conditions in the DCE experiment (4). This allows to achieve 
the same regime of velocities and accelerations, for the mir¬ 
rors as well as for the qubits, and it enables the possibility 
of joint qubit-cavity motion. Using realistic parameters, we 
can consider trajectories of constant simulated acceleration 
A = 10 15 to/s 2 during 1 ns 0, or harmonic motion with 
maximum accelerations of .4 = 10 17 m/s 2 and maximum 
velocities of around i>/4, where v = 1.2 x 10 s m/s is the 
group velocity of the electromagnetic field. 

We propose now the detection of simulated cavity- 
enhanced acceleration radiation. In this sense, a two-level 
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FIG. 3: (Color online) Expected value of a z and average photon 
number over evolution time, for a numerical simulation of the system 
dynamics (red continuous line) and for the exact anti-JC model (blue 
dashed line). We consider a resonant qubit (u; = ai q = 27r x 4 GHz), 
the qubit-cavity coupling constant is given by g(t) = go cos[7r/2 + 
7t/2 cos(wdt)], with go/hi = 0.01, and cod = 2w. The upper 
plots show unitary dynamics, while the lower ones include dissipa¬ 
tive qubit decay characterized by -y/2n = 400 KHz. 


atom in its ground state accelerated along a single-mode cav¬ 
ity has a certain probability of emitting a photon and jump to 
its excited state, even if the cavity is in the vacuum state □D- 
This is due to the fact that the acceleration activates the coun¬ 
terrotating terms of the Hamiltonian. The cavity accelera¬ 
tion radiation is characterised by the absorption and emis¬ 
sion coefficients f? a bs,em> which are related with the tran¬ 
sition amplitudes between | g 0 ) (qubit in the ground state, 
field in the vacuum) and |e 1 ), (qubit excited, one photon), 
Rem = | f^ 1 dt(el\Hi(x q , t)\g0)\ 2 and between \g2) and 

\e 1 ), R a bs = | fl 1 dt(el\'Hi(x q , i)|g 2 )| 2 , where we write Hi 
in the interaction picture. Up to second order in perturbation 
theory, we can relate these magnitudes with observables such 
as the qubit excitations (cr z ) = 1 — R em and the average num¬ 
ber of photons ( n) ph = R em ( 1 + 2 R abs + R em ). 

In the case of a uniformly accelerated atom crossing the 
cavity in the longitudinal direction, the relativistic description 
of the atomic qubit trajectory is given by 


(t) = -Jcz + AH*. 


A 


Accordingly, from Eq. ([2]) we have that 


(5) 


Rem — 


ftl 2 

dtg cos [kx q (t)\e~ iulq t e~ iut e~ lt 


( 6 ) 


where we introduced the atomic qubit spontaneous decay rate 
7 . Notice that the observation of such a phenomenon with 
current quantum optical tools is unachievable in the lab, due 


to the extreme acceleration needed in order to produce a de¬ 
tectable signal 011 - 

Let us then examine specific physical phenomena that could 
be simulated with the circuit QED architecture of Fig. [2] 
For uniform accelerations, from Eqs. 0 and <[5]>, we have 
^ = Wc 2 + _4 2 t 2 . For a resonant qubit, given the cavity 
frequency ix q = w = 27r x 4 GHz, the maximum achiev¬ 
able coupling strength and acceleration are go/io = 0.01 and 
A = 10 15 m/s 2 , respectively. Accordingly, from Eq. ([ 6 ]», the 
probability of qubit excitation due to relativistic motion is of 
the order of R em ~ 10 -4 after one trip. In order to gener¬ 
ate detectable excitations with realistic accelerations and cou¬ 
pling strengths, the qubit flying time into the cavity must then 
be increased. Hence, a natural possibility is to consider a qubit 
trajectory given by harmonic oscillations. From the Hamilto¬ 
nian of Eq. ([3J and the relation in Eq. ([4|i, we note that such 
dynamics can be simulated taking the external flux / to be os¬ 
cillating with frequency uJd and amplitude A/, around a fixed 
offset /o. 


f(t) = fo + A/cos(w d f) = kx q (t). (7) 


In this case, the absorption and emission coefficient read 


R e 


,abs — 


dtg(t) e ~ iUgt e :fiut e -' yt 


'to 


( 8 ) 


where g(t) = g 0 cos [k (x 0 + Aa; cos u>d t )]. Here, x 0 = 
f 0 /k is the equilibrium position, while Aa; = A f /k and u>d 
are the amplitude and the frequency of the oscillations, respec¬ 
tively. 

Notice that, when g(t) follows harmonic oscillations at the 
sum of the frequencies of the cavity and the qubit, the counter¬ 
rotating transitions \g 0 ) -4 |el), |e 1 ) —> |g 0 ) become res¬ 
onant, and the anti-Jaynes Cummings (anti-JC) model is ob¬ 
tained. This is the case of a qubit crossing the cavity with 
no acceleration and at the double of the speed of light in the 
material. This condition makes such dynamics not imple- 
mentable with current superconducting technology. However, 
if we set fo = 7t/2, A/ = 7r/2 and u>d = 2co, the system 
behaves as the anti-JC model with effective coupling strength 
slightly larger than go. This can be understood by developing 
the coupling time-dependence in a power series, where only 
frequency components rotating at 2 u> have non-negligible ef¬ 
fect on the system dynamics.The effective motion of the qubit, 
when coupled to the second resonant mode of a A/2 resonator, 
is given by [see Eq. (|4]i] x q {t) = L/2 [1 + cos(w,j t)]. Hence, 
such a simulation corresponds to a qubit oscillating between 
the center of the resonator and one of the two TLR capacitors. 

In Fig. [3] it is shown a comparison between numerical sim¬ 
ulations of such dynamics and the exact anti-JC model. The 
initial state for the simulation is the vacuum, i.e., the qubit in 
its ground state and no photon in the resonator. Regardless 
of the similarities between these models, we have found a pa¬ 
rameter regime showing single-excitation transfer between the 
cavity and the qubit. Interestingly, by introducing dissipation 
in the system one can force the anti-JC dynamics out of the 
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FIG. 4: (Color online) Numerical simulations for a resonant qubit 
(uj q = oj = 2n x 4 GHz) and decay rate y/2n = 400 KHz. The 
coupling strength is given by g(t) = gocos[-7r + A/ cos(wdf)], 
with go/ui = 0.01. Here, A/ = 7r (red continuous line), A/ = 
0.9-7T (green dashed line) and A/ = 0.8-7T (blue dashed-dotted line). 
The value A/ = 7r corresponds to a qubit bouncing back and forth, 
spanning all the cavity length. 


single-excitation subspace, as shown in Fig. [3] This results in 
photon growth induced by the joint effect of counterrotating 
terms and dissipation, as studied in Ref. (28l . 

On the other hand, if the frequency of the oscillations 
matches half the sum of the cavity and qubit frequencies, the 
counterrotating transitions are no longer resonant and experi¬ 
ence similar behavior as the rotating ones. Therefore, the tran¬ 
sition \e\) —> |p0) has a similar probability as |e 1 ) —> \g2 ), 
resulting in a parametric photon creation and amplification. In 
Fig. 0 we observe this phenomenon when /o = 7r, A/ = n 
and u>d = ui = oj q . Such configuration corresponds to a qubit 
bouncing back and forth between the two ends of a A/2 res¬ 
onator. This photon generation can be related to the DCE (29), 
this time related to the relativistic motion of the qubit. 

Finally, we point out that the proposed model could be 
implemented with other quantum technologies that display 
high-frequency tunability properties. For example, the sharp 
frequency dependence of giant artificial atoms ED may be 
also used to tune the effective coupling strength with negli¬ 
gible effect on the qubit frequency. Using hybrid technolo¬ 
gies ED, the difference in frequency scales could be exploited 
favourably to increase the ratio between the tuning frequency 
and the system characteristic frequencies. 

In summary, we have shown that ultrafast modulation of 
the coupling strength between a superconducting qubit and a 
single mode of a superconducting resonator mimics the effec¬ 
tive motion of the qubit at relativistic speeds. When the qubit 
follows an effective oscillatory motion, we find two different 
regimes. Depending on the oscillation frequency the system 
resembles the anti-JC dynamics or experiences unbounded 
photon generation. Our proposal provides an essential tool 
for the quantum simulation of relativistic QFT and, thus, the 
interplay between quantum physics and relativity. 
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